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Abstract 

We have employed three-dimensional numerical simulations of the impact of Comet 
Shoemaker - Levy 9 (SL9) on Jupiter and the resulting vapor plume expansion using Smoothed 
Particle Hydrodynamics (SPH) method. An icy body with a diameter of 2 km can penetrate to an 
altitude of -350 km (0 km = 1 bar) and most of the incident kinetic energy is transferred to the 
atmosphere between -100 km to -250 km. This energy is converted to potential energy of the 
resulting gas plume. The unconfined plume expands vertically and has a peak radiative power 
approximately equal to the total radiation from Jupiter's disc. The plume rises a few tens of 
atmospheric scale heights in -10 2 seconds. The rising plume reaches the altitude of -3000 km, 
however, no atmospheric gas is accelerated to the escape velocity (-60 km/s). 

INTRODUCTION 

Fragments of Comet Shoemaker-Levy-9 (SL9) are predicted to impact Jupiter in July 1994 
[Shoemaker et al., 1993; Yeomans and Chodos, 1993]. Observations using the Hubble Space 
Telescope indicate the maximum fragment diameters of -4 km, and -10 fragments are -3 km in 
diameter [Weaver et al., 1993]. The total energy released to the Jovian atmosphere upon impact of 
all of the comet fragments (CF's) will be 10 29 to 10 31 erg. 
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ENTRY AND HYDRODYNAMIC INSTABILITIES 

We have modeled the impact of CFs onto Jupiter and the resulting plume expansion using 
the Lagrangian method called Smoothed Particle Hydrodynamics (SPH) [Gingold and Monaghan, 
1977]. A Tillotson equation of state for ice was used for the cometary material [O’Keefe and 
Ahrens, 1982] and an ideal gas with y = 1.4 is employed for the atmosphere. We assume the 
atmospheric structure observed by Voyager and extend our model adiabatically into the interior of 
Jupiter. We performed calculations for two sizes of CF's, 2 and 10 km in diameter with an impact 
velocity of 60 km/s and impact angle of 40° from the zenith. Detailed description of our 
calculations is given in Takata et al. [1994], These sizes characterize the range of maximum sizes 
and energy of the Comet SL9 fragments as perceived since discovery [Chapman, 1993a; 

Chapman, 1993b; Scotti and Melosh, 1993; Sekanina et al., 1993; Weaver et al., 1993]. 

When the bolide enters Jupiter's atmosphere it induces a high temperature shock wave in 
the gas in front of the bolide. The radiation and convective flow of the shocked atmosphere around 
the front melts and ablates material from the bolide. The high speed flow of the atmosphere 
tangential to the melted material drives surface-wave Kelvin-Helmholtz (K-H) instabilities that 
grow and shed melt droplets into the flow behind the body. These particles are micron-sized, as the 
fastest growing wavelengths of the K-H instabilities are microns in length. The relative mass loss 
for large impactors is small [Bronshten, 1983]. Wavelengths much shorter -1 micron are 
suppressed by the viscosity of the melt and do not grow —this is seen in the turnover in the curves 
in Figure 1 calculated using the formulation of Keith and Banks [1990]. 

As the bolide penetrates further into the atmosphere, dynamic pressures within the shocked 
material in front of the bolide exceed the strength of the impactor. This results in fracture and then 
subsequent radical deformation of the bolide. The compression occurring during the passage of 
the fracturing wave is dependent upon R/£, where R is bolide radius and l is atmospheric scale 
height The characteristic time scale of the process of fracture wave propagation is, tf w = 2R/Us, 
(see Figure 2a). The aerodynamic stresses acting on the bolide are greatest at the center of the 
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fragment and give rise to the spreading and flattening of the bolide [Melosh, 1989; Vickery, 1993; 
Zahnle, 1992; Zahnle and MacLow, 1994]. The relative diameter was found in SPH calculations 
to increase from unity to maximum of ~4 for a 1km radius impactor and -1.7 for a 5 km radius 
impactor (see Figure 2b). However, MacLow and Zahnle [1994] point out on the basis of 
calculations conducted by assuming the bolide acts as a precompressed ideal gas, that the present 
calculations (which yield motions similar to hydrodynamic instabilities) may be spatially 
unresolved and the above scale lengths may be too large. However, this type of distortion and 
flattening has been observed in laboratory experiments of breakup of liquid droplets by gas streams 
[Ranger and Nichols, 1970]. We infer that such hydrodynamic instabilities prevent the projectile 
from expanding to the effective diameters that are as large as -20 times the undisturbed diameter 
predicted by models based only on inertial confinement [Zahnle, 1992], 

Besides K-H instabilities, Rayleigh-Taylor (R-T) instabilities occur and these are driven by 
the deceleration of the bolide and by the aerodynamic forces acting near the center of the front of 
the bolide. The particles ejected in Figure 2a are believed to result from R-T instabilities. The time 
for R-T instabilities to grow relative to the characteristic deceleration time, tdeceL is given by: 

tR-T'tdecd = {[R/(2^f)](X/R)} 172 (1) 

where 

f = in {(2 / (3C d ))(p b R / (p a ^)} (2) 

and X is the instability wavelength, Cd the drag coefficient, pa and pb the densities of the 
atmosphere and bolide. Ivanov et al. [1992] observed that upon passage of impactor through the 
Venusian atmosphere that surface instabilities grow at the front of the projectile and limited the 
distortion of the projectile. Recently, Boslough et al. [1993] observed similar phenomena in their 
calculation of SL9 impact on Jupiter. 

Instabilities are driven by velocity differences between the atmosphere and the bolide 
surface and are most prevalent near the edge of the front face of the bolide. The time for K-H 
instabilities to grow relative to the bolide characteristic deceleration time is given by 

tK-H/ldecei = (!0/(uV3))(llp/ri a ) 1,2 (R/<) 2 ' 3 (X/R)(l/f) (3) 
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where rj a and rjb are the viscosities of the atmosphere and the bolide. 

The final breakup of the bolide occurs when the instabilities having wavelengths on the 
order of the size of the bolide (X/R ~ 1) grow to large amplitudes. Referring to Figure 1, we see 
that growth rates for large scale K-H instabilities are on the order of 1 s’ * for bolides 1 km in 
radius and on the order of 10 1 s'* for 1 to 5 km radius impactors. The growth rates for R-T 
instabilities are also comparable to these values. The present SPH calculations show disintegration 
occurring at time scales predicted from these growth rates (Figure 2). Shown in Figure 2c is final 
disintegration and the stopping of the disintegrated fragments. 

We find that CF's can penetrate to an altitude of -350 km(~200 bar), in the case of an initial 
diameter, D = 2 km and to the altitude of -550 km (-800 bar), in the case of D = 10 km. The 
energy transfer from the CF to the atmosphere occurs mostly in the altitude range from -100 km 
(-10 bar) to -250 km (-100 bar) in the case of D = 2 km and -350 km (-200 bar) to -480 km (-500 
bar) in the case of D = 10 km (Fig. 3). 

THE ATMOSPHERIC PLUME 

We subsequently carried out the calculation of a plume expansion (Fig. 4). We distributed 
-10 5 atmospheric particles in model atmospheric box of 500 km width extending from 350 to -400 
km altitude in the case of the CF of D = 2 km. The energy density of deposition shown in figure 3 
is placed on atmospheric particles along the trajectory. In -10 2 seconds after the impact, the 
fireball rises a few tens of atmospheric scale heights and the unconfined plume expands vertically 
rather than horizontally in an inhomogeneous atmosphere. Buoyancy forces result in the upward 
motion of the plume. Some 25% of the energy deposited in the deep atmosphere can be 
transported above 100 km in 10 2 seconds and most of the energy is brought up above 100 km in 
-200 seconds. A mass of atmospheric gas equal to 10 times the CF mass is transported from 
below the cloud deck to above 1000 km in 10 2 seconds and a total of -20 times the initial CF mass 
is elevated above 100 km. The rising plume achieves an altitude of -3000 km, however, no 
atmospheric gas is accelerated to escape velocity. The plume transport results in the vertical mixing 
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of deep atmospheric constituents, such as NH3, H2O, and CH4, in addition to the vaporized 
cometary materials. Subsequent condensation of both deep atmospheric and cometary materials 
may be a source of dust condensates in the upper atmosphere. 

RADIATIVE SIGNATURES 

The radiation upon entry from the heated atmosphere and vaporized cometary media within 
the temporary conical cavity will be multiply Rayleigh scattered by the H 2 and He of Jupiter's 
atmosphere, and obscured by the several Jovian cloud decks. Moreover, the radiant flux may also 
be reflected from one or more Galilean satellites depending upon their position during impact. The 
observed color temperature during entry of a cometary fragment into Jupiter’s atmosphere will be 
in excess of 1()4 K. The details of the spectrum will depend on the degree of ionization which 
occurs in the bolide material and in Jupiter's atmosphere [Chevalier and Sarazin, 1994; Zahnle and 
MacLow, 1994]. In any case, the ~10^ K temperature seen in the impact flash during bolide entry 
will occur (because of radiative cooling) for only a minute or so. 

Because the optical radiation from the (~10sec) entry of a SL9 fragment into Jupiter will 
emit radiation anisotropically , we chose to model only the radiation from the second SPH 
calculation describing the plume. Moreover, the entry flash may be partially obscured in many 
directions by the asymmetry of the entry hole, opacity of ionized gas, and the NH 3 , NH 4 SH, and 
H 2 O cloud layers. 

Starting with plume calculations of 2-and 10-km SL9 fragments, we have considered each 
particle as a greybody radiator disc and calculated the total normal radiative power as a function of 
wavelength at a series of times (e.g. Figure 5). Since we are not accounting for the absorption of 
CH 4 or Rayleigh scattering from H 2 , the present calculations are an upperbound to the actual 
radiating power. However, most of the hot materials, shown in Figure 4 are at, or above, the 
elevations such that there will be minimal absorption effects of clouds. We observe that in the case 
of the 2 km comet nucleus fragment impact that its particle velocity is such that it should achieve an 
altitude of 3000 km within several minutes. 
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PREDICTIONS 

In July 1994, the Galileo spacecraft will be in a good position (1.6 AU away from Jupiter) 
to observe with the Solid State Imaging Experiment [Russell, 1992] the plumes from SL9 fragment 
collisions directly. The flux emitted by the plume associated with a 2-km diameter fragment is 
some 30% of the solar visible flux reflected by Jupiter. It is comparable to both Jupiter's reflected 
solar flux near 1 micron and an order of magnitude or more higher than Jupiter’s near-infrared 
reflected flux, and Jupiter’s average 5- and 7.8-micron thermal emission. Using our 2-and 10-km 
results, we infer that the plume from impact of a 2-km diameter fragment yields a radiant power 
equivalent to 3% of Jupiter's total flux and some 10% of Jupiter's average thermal emission at 5 
and 7.8 microns. Thus, observations by Galileo instruments are technically feasible for a large 
number of the fragments if they are as dense as we have assumed. 

Direct earth-based observations of the impact-induced plume may be possible as it rises 
above Jupiter's horizon as seen from the earth (prior to dawn at the impact point). However, 
observation of the inner Galilean satellites brightnesses as a function of time may detect a reflection 
of the plume illumination superimposed on the solar illumination. The plume generated by a 1-km 
diameter fragment would increase Io's visible flux by some 10% in the visible and by 1% or less in 
the near infrared. The plume generated by the 0.5-km fragment illuminates Io only 0.01-0.03% 
more than the sun between the visible and near infrared. However, at wavelengths of 5 microns 
and longer, the 0.5 km fragment induced plume illuminates Io some 1-3% more than the sun. 
Clearly, them most favorable spectral band for indirect detection of impart-induced plumes is in the 
near or middle infrared. 

For SL9 fragments with diameters <1 km, impact-induced plume brightnesses are expected 
to be decreased by plume and atmospheric opacity effects [Chevalier and Sarazin, 1994], In these 
cases, the wavelength of peak spectral power is likely to be diagnostic of fragment size. Our 
present study indicated the wavelength of maximum plume intensity, Xm, and fragment diameter, 
D, are related by the Wein's law type relation: X m D3 ~ 5 x 10^ cm'*. 
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Fig 1 Growth rate of Kelvin-Helmholtz instabilities as a function of the diameter of the impactor 
divided by the instability wave length. The curves are for impactor diameters of 1, 2, and 10 km 
and the depth into the atmosphere is 300 km relative to the 1 bar level of Jupiter. 

Fig. 2. Smoothed-Particle -Hydrodynamic (SPH) calculations of the penetration of a fragment of 
Comet Shoemaker-Levy-9 fragment into Jupiter's atmosphere. The diameter of the bolide is 10 
km and the angle of entry is 40° from the zenith. Shown are the density fields at various times 
relative to the time the bolide penetrated to the 1 bar altitude. Dots indicate planetary atmosphere 
and triangles bolide SPH particles. The times are 2.0 , 5.0, and 21 sec., respectively . 

Fig. 3. The fractional energy deposition of the initial kinetic energy of the comet as a function 
altitude in the case of 2 and 10 km diameter-cometary fragments. The results of SPH calculation 
are shown by the solid lines, and the results of meteoric ablation model with the ablation 
coefficients, o, of 10 -9 and 10- 8 s 2 /m 2 [Bronshten, 1983] for D = 10 km, are shown by dashed 
lines. 

Fig. 4. The pressure of atmospheric particles upon plume expansion of D = 2 km cometary 
fragment after 75.2 sec (after Ahrens et al. [1994]). 

Fig. 5 Calculated spectra from Shoemaker-Levy-9 impact-induced plumes for 2 and 10 km 
diameter fragments. The peak power from the 2 km fragments radiates at the approximately same 
power level as the total Jupiter disk. The total power from the sun, considered as a point source is 
shown, for comparison. 
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